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INTRODUCTION 

Numerous  applications  arise  where  it  is  desirable  to  apply  the 
spectral  dispersive  properties  of  holographic  elements  in  order  to  perform 
a spectral  decomposition.  To  design  holograms  for  this  type  of  application 
in  most  cases  does  not  require  an  in-depth  understanding  of  the  holographic 
process.  The  following  is  a discussion  of  the  basic  concepts  and  techniques 
of  sufficient  detail  to  permit  the  reader  to  design  and  use  holographic 
optics  for  a number  of  spectral  analysis  applications.  Included  is  a 
calculator  ray  trace  program  for  use  on  the  Hewlett  Packard  9820A  calcu- 
lator. This  program  can  serve  as  a design  aid  for  a large  class  of 
applications. 
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SECTION  II 


BASIC  VECTOR  RELATIONSHIPS 

The  field  emanating  from  a point  source  can  be  represented  by. 


Aelk*  r (1) 

where  A is  the  amplitude  of  the  field,  k*  - (2ir/X)k  is  the  propagation 
vector,  and  r1 is  the  vector  from  the  point  source  to  the  point  of  obser- 
vation. The  phase  relationship  is  described  by  the  vector  dot  product 


k • r (2) 

In  the  construction  of  a hologram  two  such  fields  kR  and  kQ  (refer- 
ence and  object  fields)  exist  simultaneously  within  the  recording  material. 
The  resultant  field  is  characterized  by  the  vector  relation 


- v • *- 

k • r o) 

where  K is  the  vector  describing  the  orientation  of  the  resultant  inter- 
ference fringes. 

The  recording  and  processing  of  the  imposed  fringe  pattern  within 
the  recording  material  results  in  a fringe  pattern  that  is  proportional 
to  this  incident  interference  pattern.  (For  a more  complete  description 
of  this  process  in  silver  halide  emulsions,  see  reference  1)  i.e., 
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where  Ky  is  Che  vector  describing  the  orientation  of  the  recorded  fringes 
and  a is  a constant  reflecting  the  exposure  and  recording  characteristics. 

The  process  of  reconstruction  consists  of  illuminating  the  processed 
hologram  (in  the  same  angular  orientation  as  recorded)  with  a wavefront 
that  duplicates  the  reference  wave.  The  resultant  diffracted  wavefront 
has  a propagation  vector  given  by, 

k - ty  + nK  (5) 

where  n is  the  diffraction  order.  If  the  wavelength  and  geometry  is 
identical  to  those  of  the  construction  process,  the  propagation  vector  will 
reduce  to, 

k - nkQ  , (l+n)kg  - nk^  (6) 

Only  the  first  term  is  of  interest  for  this  usage.  The  first  diffracted 
order  then  results  in  a wavefront  that  is  identical  to  the  original  object 

beam. 

If  a wavelength  variation  is  included  between  the  reconstruction  and 
construction,  the  vector  relationship  becomes 

k(2r/Xx)  - (2w/XI)ltR  + n(2ir/XE)  K 

k ■ kg  + n(Xj/XE)  K.  (7) 

where  * indicates  the  unit  vector  and  the  subscripts  E and  I refer  to 
exposure  and  illumination  conditions.  Thus  for  a given  angular  direction 
of  the  reconstructing  field,  the  angle  of  the  reconstructed  ray  is 
dependent  on  the  illumination  wavelength. 
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This  analysis  considers  the  fringe  spacing  to  be  a constant  over 
the  area  of  reconstruction.  Through  superposition,  this  can  be  extended 
to  more  complex  holographic  structures  and  Is,  In  fact,  quite  general. 

The  important  result  is  that  the  change  in  propagation  direction  of  the 
illumination  that  is  induced  by  the  hologram  is  proportional  to  the  wave- 
length of  the  radiation.  This  result  is  illustrated  in  figure  1.  This 
figure  shows  the  change  in  focal  position  that  results  from  various  wave- 
lengths of  illumination  when  the  holographic  grating  is  placed  in  front 
of  a simple  lens.  From  this  figure  it  is  obvious  that  the  plane  containing 
the  on  axis  focal  points  for  the  various  wavelengths  is  curved.  This  will 
result  in  a design  limitation  if  image  resolution  is  to  be  maintained  over 
a significant  range  of  wavelengths. 

One  other  feature  that  must  be  taken  into  account  in  the  design  is 
the  fact  that  the  multiple  diffraction  orders  may  overlap  either  due  to 
spectral  spread  or  spatial  extent.  The  overall  range  of  wavelengths  is 
limited  only  by.  the  angular  range  of  the  optics  and  the  transmission 
characteristics  of  both  the  optics  and  the  holographic  materials.  A 
considerable  amount  of  flexibility  exists  due  to  the  ability  of  varying 
the  spatial  frequency  of  the  holographic  fringes  in  the  construction 
process. 
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Hologram  Curved  Focal  Plane 
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SECTION  III 
CONCLUSION 

The  proceeding  discussion  is  meant  to  be  only  a first  look  at  the 
use  of  holographic  gratings  as  a tool  to  perform  spectral  decomposition. 

The  discussion  gives  some  of  the  basic  vector  relations  that  are  used  in 
the  calculator  ray  trace  program  of  the  attached  Appendix.  This  discussion 
together  with  the  ray  trace  program  will  allow  the  reader  to  perform 
sufficient  design  analysis  to  answer  the  first  questions  associated  with 
the  use  of  holographic  elements.  For  a more  detailed  treatment  of  the 
theory  of  holographic  optics,  the  reader  is  referred  to  any  of  the  current 
texts  on  the  subject  or  the  three  dissertations  listed  as  references. 


APPENDIX 


CALCULATOR  RAY  TRACE  PROGRAM  FOR  SPHERICAL 
DIFFRACTIVE,  REFRACTIVE  AND  REFLECTIVE 
OPTICAL  ELEMENTS 
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Introduction 


Over  the  past  several  years  programs  have  been 
written  to  simulate  various  aspects  of  hologram  perf- 
ormance. These  programs  vary  extensively  In  the  degree 
of  complexity  depending  on  the  flexibility  and  thorough- 
ness required.  Several  programs  for  use  with  the  table 
top  (Hewlett-Packard  9820  ) calculator  have  existed  in 
DHO  alone.  These  programs  are  typically  written  for  a 
particular  problem  and  for  this  reason  they  are  not 
general  enough  for  widespread  usage.  This  is  an  attempt 
to  write  and  document  a relatively  simple  program  that 
is  sufficient  to  be  used  for  the  analysis  of  a large 
number  of  optical  problems.  The  program  will  handle 
refractive  elements  of  spherical  shape  as  well  as 
reflective  spherical  elements  and  holographic  elements. 

The  program  will  provide  a ray  trace  analysis  of  an 
optical  system  that  is  a combination  of  refractive  or 
reflective  elements  together  with  one  holographic  lens. 

The  program  uses  Snell's  law  to  calculate  the  ray 
direction  across  the  refractive  Interfaces,  and  the 
vector  relationship  to  determine  the  ray  direction 
through  diffractive  elements.  For  reflective  elements 

8 
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the  index  of  refraction  and  the  surface  to  surface 
distance  must  undergo  a sign  change  at  the  reflective 
Interface.  This  will  be  illustrated  by  an  example  in 
the  ensuing  discussion. 

The  following  discussion  gives  a description  of 
the  program,  Its  Input  requirements,  the  memory  alloc- 
ations, the  calculations  performed,  examples  of  each 
type  of  element,  and  a printout  of  the  program. 

The  calculations  are  shown  In  mathematical  form  together 
with  their  equivalent  statement  In  calculator  form. 

This  is  provided  for  ease  In  trouble  shooting  as  well 
as  program  modification.  The  program  Is  written  for 
use  on  the  9820A  Hewlett-Packard  calculator  In  conjunc- 
tion with  the  9862A  Plotter,  9866A  Printer,  and  the 
9865A  Cassette  Memory. 
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INPUT  PARAMETERS 


SCALEJ  Xmin’  Xmax’  W Ymax  * 


NON-HOLOGRAM  SURFACES: 


Radius  of  curvature  of  the  surface. 

( Positive  If  center  of  curvature  Is 
to  the  right  of  surface. )(  0 « flat  surf.) 


On  axis  distance  to  the  next  surface. 


Index  of  refraction  of  the  medium 
following  the  surface. 

Clear  Aperture.  ( Radius  of  the  active 
area  of  the  surface.) 


INPUT  PARAMETERS 
(continued) 


HOLOGRAM  PARAMETERS: 


R/D  Diverging  wave  focal  distance  from  the 
hologram  surface. 


“/D  Diverging  wave  angle  with  axis. 

(Positive  If  point  source  Is  above  the 
axis) 


R/C  Converging  wave  focal  distance  from  the 
hologram  surface. 


«/ C Converging  wave  angle  with  axis. 

(Positive  If  point  of  convergence  Is 
below  axis) 


a/I  Illumination  wavelength. 


a/E  Exposure  wavelength. 
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INPUT  PARAMETERS 
(continued) 


INPUT  RAY  BUNDLE  PARAMETERS: 


CENTRAL-Yj  Y position  of  the  central  ray 
In  the  Input  bundle. 


DELTA-Y.  Spacing  between  the  parallel 

1 Input  rays. 


A12  Angle  of  rays  between  the  first 

two  surfaces.  (Positive  If  the 
ray  has  a negative  slope.) 


MEMORY  ALLOCATION 


REGISTER 

A 

B 

C 

X 

Y 

Z 

R ( 0 ) 
R(l) 
R(2) 
R(3) 
R(4) 
R(5) 
R(6) 
R(7) 
R(8) 
R(9) 


INFORMATION  STORED 

Surface  number  of  hologram. 
(If  none,  A«0) 

xm1n 

Xmax 

Ym1n 

Ymax 

Index  Register  for  Surfaces 
Hologram  R/D. 

Hologram  «/D. 

Hologram  R/C. 

Hologram  «/C. 

Hologram  LAMBDA/ I. 

Hologram  LAMBDA/E. 

Central  Yj . 

Delta  Yj  (From  R(498)). 
Interim  storage  used  In  cal 
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MEMORY  ALLOCATION 
continued 


REGISTER  INFORMATION  STORED 

R(IOZ)  Radius  of  curvature  of  surface  Z. 

R( 1 OZ+1 ) On  axis  distance  to  surface  Z+l . 

R(10Z+2)  Index  of  refraction  between  Z and 

Z+l . 

R ( 1 OZ+3)  Radius  of  clear  aperture  of  Z. 

R(1 OZ+4)  Angle  of  Incoming  ray  at  plane  Z. 

R(10Z+5)  X position  of  Z plane  on  axis. 

R(10Z+6)  X distance  between  Intercepts  of 

planes  Z-l  and  Z. 

R( 1 0Z+7 ) X value  at  Intercept  point  on 

plane  Z. 

R(1 OZ+8)  Y value  at  Intercept  point  on 

plane  Z. 

R(10Z+9)  Angle  of  normal  to  surface  at 

Intercept  point. 

R(498)  Delta  Y1 

R{ 507 ) Interin  calculations  In  Holo.  sub. 

R(508)  " " " " 

R ( 509 ) " " “ 


DEFINITION  OF  TERMS 


XZQ  On  axis  x component  of  Zth  surface. 

d7  7.,  X distance  between  the  Intercept  on 

* the  Z and  Z+l  surfaces. 

*Z  Z+l  ^ngl e with  respect  to  the  x axis  of  the 

* ray  between  surfaces  Z and  Z+1 . 

"RZ  Angle  of  the  normal  to  the  Z^  surface 

at  the  point  of  Intersection. 

X7  x component  of  the  Intersection  point 
on  the  Z^h  surface. 

Yz  y component  of  the  Intersection  point 
on  the  Zth  surface. 

n7  Index  of  refraction  of  the  medium 
following  the  Zth  surface. 


x.1.  x component  of  the  center  of  curvature 
L of  the  Zth  surface. 


Calculations  for  a flat  surface: 


i 

| 


Figure  A-l.  Flat  Surface  Geometry 

™t  + X0  ’ X20 

R<Z+1)  ♦ R(Z+5)  -*►  R(Z+15) 

X20  ’ X1  " d01 

R(Z+1 5)  - R(Z+7)  R(Z+1 6) 

IF  NEXT  SURFACE  IS  FLAT 
IF  R(Z+10)  ■ 0 

Y1  “ d01  t#n*12  " V2 

R(Z+8)  - R(Z+1 6)tanR(Z+4)  R(Z+18) 
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0 -*•  R( 2+1 9) 

Sin-1  ({n^ngjSIn  « 12)  * “23 

Sln-^RU+aJCSInRtZMjJ/RtZ+lZ))  -*-R(Z+14) 

GO  TO  NEXT  SURFACE 
JMP  8 


l 

i 

\ 

f Calculations  for  a spherical  surface: 

The  calculations  for  the  ray  intercept  with  a 

. spherical  surface  require  the  simultaneous  solution 

( 

t of  the  equation  describing  the  ray  path  and  that  of 
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Figure  A-2.  Geometrical  Layout  for  the  Spherical 
Surface  Calculation. 


the  spherical  surface.  The  equation  of  the  ray  Is, 


where 


y * mx  + b 


m * -tan 


Then 


b * Y'j  + Xjtan  «j2' 


-x  tan  + (Yi  + X^tan 


The  equation  for  a circular  Intercept  whose  center 
Is  located  on  the  x axis  Is  of  the  form; 


where 


( x - x;,)2  + y2  * Rno2 


x12  * x20  + R02 


Is  the  center.  Thus 


y2  * Rq22  - x2  + 2xj'2x  - x.' 
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Setting  these  equations  for  y equal  and  solving  for 
the  Intercept  point  X2  and  letting 

2 

B * (xj2  + (T^+  X^tan  “i2^Cos  “l2 

and  « a o 

C * *x12  " RdI  + *Y1+  Xltan  *12'  )Cos  “12 

we  have  the  familiar  equation 

2 

x2  - 2Bx2  + C « 0 


Thus 


X2  » B t ( B2  - C )1/2 


and 

Y2  * Yi  + (xi“  X2)tan  -12. 

The  sign  Is  determined  by  ; 

rD2  ^ ^02 

l«D?l  " ’ 

This  solution  Is  implemented  in  the  program  by 
the  following  steps. 

X20  + RD2  + X12 
R(Z+15)  + R(Z+10)  -*•  R9 

Calculation  of  B 


E ! 


CoszR(Z*4)(R9  » (R(Z+8)+R(Z+7)tanR(Z.4))tanR(Z+4))  -»R8 


Calculation  of  C 

Cos2R(Z+4)(R92-R(Z+10)2+(R(Z+8)+R(Z+7)tanR(Z+4))2>-*R7 


B - R02((B2-  C)/Rd22)  1/2  - x2 
R8  - R(Z+10)( { R8 2 - R7)/R(Z+10)2)1/2-*.  R(Z+17) 
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Y1  - (X2-  X1  )tan  «12  * Y2 

R( Z+8)  - (R(Z+1 7) - R(Z+7))tanR(Z+4)  -►  H( Z+l 8 ) 

stn-'(  r2/RD2)  - -.R2 

S1n-1(R(Z+18)/R(Z+10))  -^R(Z+19) 

“R2+  Sin-1  ((n1/n2)(Sin(  “12-  *23 

R(Z+19)+Sin-1((R(Z+2)/R(Z+12))(Sin(R( Z+4) -R(Z+19))/R(Z+12))) 
-*•  R(Z+14) 

GO  TO  NEXT  SURFACE. 
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Calculations  for  a holographic  element: 


i 


! 


The  calculations  of  the  ray  path  through  a 
holographic  element  are  based  on  the  following 
vector  relations.  The  fringe  pattern  recorded  in 
the  hologram  will  have  a localized  fringe  distrib- 
ution whose  orientation  is  described  by  the  vector 
K where, 

K » - 72  . 

kj  and  k2  are  the  wave  propagation  vector  for  the 
two  fields  used  in  the  construction  of  the  holo- 
gram. 

The  vector  K can  then  be  used  to  determine  the 
direction  of  the  diffracted  ray  upon  reconstruction. 

iR  ' ri  • 7 

where  kR  and  kj  refer  to  the  reconstructed  and 
Illumination  rays  respectively.  The  program  uses 
the  direction  sines  for  the  wave  vectors  for  the 
calculations . 

The  Implementation  is  as  follows; 

Direction  sine  of  the  diverging  field  Is 

S1n0l 2D*  n/n1)s1r»(tan-l((Y1+R/DS1n«/D)/R/DCos./D)) 


j 

j 

J 

3 

! 


J 

5 


i.l 
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Sin  tan-1((R(Z+8)+R0S1nRl)/R0CosRl)/R(Z+2)  R508 

The  direction  sine  of  the  converging  field  is 

Sinfl12C=  (1/n-j  )S1ntan_1  ((Y^+R/CSin<*/C)/R/CCos<*/C) . 

Sin  tan"1 ((R(Z+8)+R2S1nR3)/R2CosR3)/R(Z+2)  ->■  R507 
Then,  the  diffracted  ray  angle  is 

9'  * Sin"1  (S1n“  j 2 / A£)(S1n0-|2D~^^n®12C^ 

Sin'1(SinR(Z+4)-(R4/R5)(R508-R507))  -*>  R509 

“12  R510 
R(Z+4)  — » R 5 1 0 

^TH1  + X10  * X20^  X20  ' X1  * d01) 

(R( Z+l ) + R ( Z+5 ) — > R(Z+1 5) ) - R(Z+7)  — R(Z+1 6) 

Y1  - d01tan9'  - Y? 

R( Z+8)  - R(Z+16)tanR509  R(Z+18) 

X1  + d01  * X2 


R( Z+7)  + R(Z+1 6) 


R(Z+17) 


Data  Output: 


The  program  Is  designed  primarily  to  provide 
a ray  trace  plot  as  the  output.  This  plot  consists 
of  a surface  to  surface  plot  of  each  of  the  rays 
In  the  prescribed  ray  bundle.  An  optional  feature 
that  exists  Is  the  plotting  of  the  surfaces  them- 
selves. This  subroutine  Is  executed  by  the 
answer  " 1 RUN  PROGRAM"  when  the  display  "SURF.?" 
is  shown.  The  plot  of  a given  ray  will  be  stopped 
If  the  Intersection  with  the  next  surface  exceeds 
the  clear  aperture  range. 

A second  form  of  output  Is  the  tabular  listing 
of  the  Input  surface  parameters.  This  output  Is 
Is  Included  for  record  keeping  and  as  a check  of 
the  accuracy  of  the  Input.  The  format  of  the 
data  printout  Is  shown  In  Figure  A-3. 
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Figure  A-4 . Sample  plot  for  a simple  lens  with  an 
Input  ray  bundle  with  spacing  of  .1 
and  an  angle  of  20°.  The  surface  param 
eters  are  shown  In  Figure  A-5. 
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Figure  A-e.  Sample  Ray  Trace  through  a Holographic 
Lens.  Input  Rays  are  the  same  as 

! Figure  A- 3. 
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3> 

I 

< 

9> 

U 

3> 


i 

l 


0: 

FXB  6 5SFG  145 
ENT  "SCALE, X MIN 
"«Bi "X  MAX" ^ Ci "Y 
MIN" i Mi "Y  MAX"  » 
YF 
1 : 

ENT  HuLO.AT  SUR 
F. NO?" i AS  IF  A = 8 i 
JMP  3F 
2: 

ENT  ” R/D  = ?" > R0i " 
ALPHA /D*? " i R1 1 "R 
/C  = ?"  i R2 » "ALPHA-' 
C=?" i R3F 

O * 

ENT  "LA  M D A / I = ? “ i 
R4 i "LAMDA/E=?" i R 
5H 
4: 

10  + ZF 
5: 

ENT  "RDs?"i RZi “T 
H ~ ? " i P <2+1.1  i “N»? 

" i R i 2 + 2 ) i " CA  = ? " i 

R ( 2+3 ) h 

b*. 

IF  R(Z+1)*05Z+10 
+2? JMP  -IF 

7 ■ 

ENT  "CENTRAL  Yl= 
?" i R?  > " DELTR-Y 1 = 
?"  i R498 , "A 12  = ?“ i 
R14F  ...... 

y • 

R 7 - ('•'  1 3 + R 1 8 i R498  + 

R8F 

9: 

0+R 1 5 + R 1 6+R 1 7 + R 1 

9 ? 1 

10: 

(R(Z+1 ) +R ( 2+5) *R 
(2+15) ) -R ( 2 + 7 ) + R 
( Z + 1 6 ) i IF  R ( Z + 1 0 
) *0  5 JMP  3F 

1 1 : 

R ( 2 + 8 ) -R (Z+l 6) 
TAN  R (Z  + 4) *R (Z+l 
8)  ? R ( Z + 7 ) + R ( 2 + 1 6 
J +R  ( Z+ 1 7 ) F 
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12: 

0+R (Z+l 9 i i ASN  (R 
(Z  + 2)  SIN  R i 2 + 4.1  / 
R • 2 + 1 2 ) J+Ri.  <£+14) 
i JMP  8F 
1 3 : 

R ( Z + 1 5 I + R f Z + 1 A i * 
R (2+1 9 ) F 

14: 

(R  ( Z + 1 9 ) + (R  (Z  + 8'J 
+R (2+7) TAN  R ( Z + 4 
) ) TAN  R (Z+4) ) 

COS  R (Z+4) t2+R (2 

+ 1 8 ) F 

15: 

(R(Z+19)  t-2-R  (2  + 
0) t2+ (R (Z+S) +R ( 

+ 7 ) T A N R ( Z + 4 ) ) ? 

) COS  R (Z+4) +2+R 
Z+l  7)  F 
16: 

R C2+3 ? ) - R (Z+l O j 
( (R ( Z + 1 8 ) 1 2 - R (Z 
17)  j /R  (2+10)  1-2) 

R ( Z + 1 7 ) F 
17: 

R (2+8) - (R (2+17) 
R( Z + 7)) TAN  R ( 2 + 

) + R ( 2 + 1 8 ) F 
1 8 : 

ASN  <R(Z+18)/R( 

+ 10) ) + R ( Z + 1 9 ) F 

19: 

R i 2+ 1 9 ) + ASN  ( R ( 
+2) SIN  (R ( Z+4) - 
( Z+  1 9 ) ) /R  ( Z + J i: ) 

P ( Z + 1 4 ) F 
20: 

IF  R (Z+l ) *05 2+10 
+ Z > IF  2=10 A 3 G S 8 
"HOLD " F 


■ x>  r-j  r j & i •+  + '-i  f'.<  r j ►- 


r 


r 


i 

21 : 

IF  RIZ+1)*0S JMP 
-111- 

22: 

SCL  BjCjXjYI- 
23: 

10-^21- 
24: 

IF  fi B S R 1 2 + 8 ) > R ( 
2 + 3)5  JMP  SI- 
25: 

PLT  R(2+7)  » F:  (2  + 3 
) P 
26 : 

IF  R(Z+1)*05Z+10 
+2 5 JMP  -2P 
27: 

PEN  5R18+R8+R185 
IF  R18£ R7+R135 
GTO  9P 

005 

ENT  "SURF. ? " j R 1 0 
00  5 IF  R 1000= 15 
GSB  "SURF . " P 
29: 

ENT  "PRINT  OUT  ? 
" > R 1 00 1 5 IF  R1001 
- 1 5 GSB  "PRINTER" 
P 

30: 

GTO  7P 


"SURF.  "P 
32 : 

28->ZP 

33: 

JMP  15  IF  ft* 0 5 IF 
2=100  + 20  5 JMP  6P 
34  s 

IF  R Z = 0 5 PLT  R(Z  + 
5) » -R (2  + 3) 5 PLT  R 
(2+5) > R (2+3) 5 
PEN  5 IF  R i (2  + 10  + 
2 ) -9  )'*0  5 JMP  -IP 
35 : 

-R (2+3) + Ri55 IF  R 
(2-9) =05 JMP  4P 

O fl  j 

PLT  R (2+5) +R2-R2 
r ( 1 - R 1 5 1 2 / R 2 1 £ ) , 
R15P 
37: 

R 1 5 + R ( 2 + 3 ) / 2 0 + R 1 
55  IF  RISiR (2  + 3)5 
JMP  -IP 
38* 

PEN  5 2+  1 0 + 2 5 IF  R 
(2-9) *05 JMP  -5P 
39: 

RET  P 


I 
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40 : 

"HOLD  “ P 
41  : 

SIN  fiTH  ((RiZ+y) 
+R0SIN  R1)/R0 
COS  R 1 j / R ( 2 + 2 ) + R 
588P 
42: 

SIN  fiTH  UR(Z  + S) 
+ R2SIH  R 3 ) / R 2 
C 0 S R 3 ) R ( 2 + 2 ) + R 

5S7(- 

43* 

flSN  (SIN  R i Z + 4 j - 
R4 (R508-R5O?) /R5 
.1  + R 5 O 9 ? R i.  Z + 4 .i  9 R 5 
1 @P 
44: 

i.  R ( Z + 1 j + R ( Z + 5 ) 4 R 
(2  + 15) ) - R ( Z + 7 ) 4 R 
( 2 + 1 6 ) P 
45: 

R ( 2 + 8 ) -R ( 2+ 1 0 J 
TflN  R509+R (Z+18) 

! R i.  Z + 7 J + R ( Z + 1 6 ) •'* 
R (2+1? ) I- 
46 : 

IF  flBS  ( R ( Z + 2 ) 
SIN  R 5 O 9 R ( Z + 1 2 ) 
->R  (Z+14)  ) i 1 5 flSN 
R ( 2+ 1 4 ) +R (2+14) P 
47: 

Z+10  + Z5GT0  101- 


" PR  INTER  "P 
49 : 

FMT  2/ 5 l.iJRT  SP 
50 : 

IP  fl  = 0JJI1P  7P 
51 : 

F M T " H 0 L 0 . H 7 " ? 
UR 7 81- 


FMT  "SURF. HO . “ , 3 
"R---D"  j 3H>  "RLPH 
H/D" , 3X» "R/C" . 3H 
* "RLPH H ■ ' C " 1 :•! X < 7 5 
I.1.IR  T 8P 

C*  . . 


FMT  " L H M J) fi / 1 " , 3X 

> “lamiih.-e  " ; mrt  s 

p 

54: 

FMT  FXD  4.0.FXI1 
11.2,  FXI:  8.2,2 
F X D 8 . 2 f 2 F X Ii  10, 
4>FKB  11.4,FXD  1 
G.2?  l.t.l  R T Sh 

c-  cr  . 

J J » 

WRT  8, fl, R 0 < R1 , R2 
j R3, R4, PSP 
56 : 

FMT  2/5  UR T si- 
57: 

FMT  'SURF . HO , " > 9 
,x, , " R Ii " ? 1 6 X . " 7H  " ? 

1 6 X j " N ’ j 1 7 X ■ " C:9  " 

> i ? W R T 81- 

17  . 

•u.1  O » 

10  + ZP 
59: 

F M T F X D 4 .0-4 
F X D 18.6?  U R T 3*2 
1 O , RZ  j R (2+  1 1 - p <• 
2+2.1  ? R (Z  + 3)  P 
66 : 

IF  r ( z + 1 1 * 0 ; 7 + 1 0 
+Z-JHP  -IP 
61  : 

FMT  4 URT  y ; 

RF.T  y~ 

62 : 

STP  P 
63: 

END  P 
R 1 1 9 9 
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